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ABSTRACT 



In Gallium detectors for solar neutrinos, the capture rate due to Gamow- Teller tran- 
sitions to excited states of ^^Ge beyond the neutron emission threshold is usually 
neglected. We make a model calculation to estimate its effect and find that this yields 
an additional contribution which may be as much as 0.4 SNU, even larger than that 
from the Isobaric Analog State in ^^Ge reached by Fermi transitions, which is normally 
included in the standard predictions. 
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Research in the area of the solar neutrino deficit problem has reached a very exciting 
stage. On one side, several groups are working on making both the theoretical as well 
as the observed neutrino capture rates for the different detectors more precise. On the 
other side, explanations of this shortfall invoking non-zero neutrino mass and mixing 
in the two- and three-flavour scenarios are being pursued vigourously. Out of the four 
detectors which have already published data, measurements using Gallium (by the 
GALLEX and SAGE groups) have the lowest threshold (0.233 MeV) and are able to 
capture neutrinos produced through all the reactions in the sun, including the most 
abundant pp neutrinos. The theoretical prediction for the total neutrino capture rate 
in the Ga detectors is 132lg SNU (Solar Neutrino Unit) using the Standard Solar 
Model (SSM) of Bahcall and Ulrich whereas it is 137i? in the SSM of Bahcall and 
Pinsonneault |0] which includes He and heavy element diffusion. To calculate the total 
capture rates one integrates the solar neutrino flux times the capture cross-section over 
the range of neutrino energy and performs summation over all the neutrino sources. 
Evaluation of the neutrino capture cross-section as a function of the energy needs 
the knowledge of the Gamow- Teller (GT) strengths to excited states of the daughter 
nucleus ^^Ge. The theoretical calculations of these GT strengths depend sensitively on 
the choice of the interaction Hamiltonian §]. The SSM predictions of ref. [0] and 
use the GT transition strengths deduced from the forward angle {p, n) reaction data 0] 
from the analogous hadronic process '^^Ga{p,ny^Ge. The neutron emission threshold 
in ^^Ge is 7.4 MeV, and, in the prevalent practice, one does not consider any v capture 
which excites states above this energy through the GT operator. This is based on the 
assumption that for such states the partial width for neutron emission is much larger 
than the 7-decay width. In this letter we go beyond this assumption and probe the 
effect that an inclusion of the above states might have on the final solar neutrino 
capture rate. To set this in perspective, the Isobaric Analogue State (IAS) in ''^Ge is 
observed at an energy of 8.932 MeV and Champagne et al. measured the upper 
limit of the ratio of the gamma decay width to the total width for the IAS to be about 
10%. Based on this, Bahcall and Ulrich [|^ add 10% of the capture rate coming from 
the Fermi strength and this turns out to yield 0.2 SNU for the solar neutrinos. But 
the n-emission from the IAS to low-lying states of ^Ge with T = 3 is suppressed by 
isospin conservation. On the other hand, for the states below the IAS {i.e., below 8.932 
MeV of excitation in ''^Ge) with T = |, connected by the GT transition, there is no 
such inhibition due to isospin conservation and hence the neutron emission probability 
is expected to be much larger. In this work we consider the contribution, albeit small, 
coming from GT-transitions to states above 7.4 MeV using information from beta- 
delayed particle emission experiments. We make a model calculation for this and to 
the best of our knowledge no earlier estimate of this possible additional contribution 
has been made. Efforts are on for the experimental measurement of the ratio of the 
gamma decay width to the total width for these states. 

The total capture rate R for ^^Ga is given hj R = T,Ri, where Ri is the capture 
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rate for neutrinos from the i^^ source. Ri is obtained by convoluting the calculated 
cross-section with (puii^l), the flux of neutrinos of energy q from the source i 0. The 
normalisation of the latter is fixed by the SSM under consideration and one has: 

'I<PuAq) (^siq) dq 

where a and as denote the absorption cross-sections as calculated from this work and 
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as given in ref. respectively. Ci is the prediction in SNU for neutrinos from the 
source according to the SSMs of ref. |l|] or p| depending on which model is used as 
the reference point. 

The capture rates a{Ei) of neutrinos from the ground state of ^^Ga to a state of 
^^Ge at an excitation energy Ei has two parts - one coming from the Fermi transition 
and the other from GT transitions - and can be written as, 

a{Ei) = apiE,) + acAEi) (2) 

where ap^Ei) and acT^Ei) are given by, 

apiE,) = 0.1{7fc'n^)-' gv'BF{E,)p,eF{Z,e) (3) 

aaAE,) = {7fc^h^)-'g\BGT{E^)peeF{Z,e) (4) 

In eqs. and gy, Qa are the vector and axial vector coupling constants, Pe and 
e are the momentum and energy of the emitted electron and F{Z, e) is the /3-decay 
Fermi function for which we use the analytic form given by Schenter and Vogel 0. 
Bp{Ei) and BcriEi) are the squared Fermi and GT matrix elements for '^^Ge states 
at energy Ei respectively. The total Fermi strength of (A^ — Z) is concentrated at 
EjAS (~ ^{Ei — EjAs)) and the net contribution is proportional to Bp{IAS). For the 
GT strength distribution the usual practice is to use BGT{Ei)s from {p,n) reactions. 
Here we are interested in the GT strength beyond the neutron emission threshold for 
which no [p, n) study results are available and we use a simple theoretical model for 
the strength distribution based on the configuration structures of the initial and final 
states which get connected by the one body GT transition operator. The model for 
the GT strength distribution is based on a formalism constructed for calculating the 
beta decay /electron capture rates in stellar conditions in ^ and used recently in 
for calculating the neutrino capture rates for the ''^Ga detectors. In this method, one 
uses pure shell model configurations to describe the ground state of Ga as well as the 
ground and excited states of ^^Ge. For example the ground state of ^^Ga is taken as 
(l/i)^ (2p3)^ for the protons and a completely filled fp shell for the neutrons. Out 
of the eight possible GT transitions: (1) 2pi{n) — * 2p3(p), (2) 2pi{n) —> 2pi{p), (3) 
l/5(n) ^ (4) l/|(n) ^ lf.{p), {5)2pi{n) -^'2ps{p), (6/2^3(71) ^'2p^{p), 

(7) 1/7 (n) 1/7 (p), (8) 1/7 (n) connecting the ground state of "^^Ga to 
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eight pure configurations in ''^Ge, (4) and (7) are completely blocked. Each remaining 
configuration has many individual states in it and the strength distribution to each 
configuration gives rise to a resonance taken to be of the Gaussian form in this model. 
The centroid energy of the strength Gaussians is given by, 

Ecent — ^Egp + /\E2—body ~\~ ^Epairing (5) 

where ^Eg.p. and AEpairing are the differences in the total single particle energies and 
pairing energies, respectively, of the excited and the ground state configurations of 
^^Ge. AE2-body is the contribution coming from the 2-body potential and is treated 
in a simplified TDA model as in [Q. The total GT strength to each configuration is 
given by the sum rule expression of the product of single particle transition strength, 
the number of neutron particles and the fraction of proton holes 0. The widths {a) 
of the strength Gaussians to different configurations are all taken equal and is a free 
parameter of the model. One finds that a = 1.55 MeV gives the total capture rate as 
132.83 SNU (137.17 SNU) close to the standard values according to ref. ^ (0). So, 
for the purpose of comparison, we use the strength width of 1.55 MeV. 

When one includes the contribution from all states beyond the neutron threshold 
of 7.4 MeV in "^^Ge the capture rate becomes 133.56 SNU (138.005 SNU), an increase 
of 0.73 SNU (0.84 SNU) over the usual SSM prediction of ref. (H). But this is 
clearly an overestimation as these states mostly decay by particle emission and not by 
gamma emission. To get a realistic estimate we use information collected from beta 
delayed particle emission experiments. With F^, being the partial widths of gamma 
and particle decay respectively, the total width F = F^ + F^. Denoting the particle 
separation energy by B^, "F^/F is zero below and remains small immediately above 
that, then increases in the space of a few hundred keV to nearly its maximum value. 



rises more slowly and finally decreases again at much higher excitations" [T^. This 
indicates that F^/F first falls slowly, then very rapidly and then slowly goes towards 
zero. We model this by a Gaussian fall-off and include that as a multiplicative factor in 
the GT strength distribution beyond 7.4 MeV. The width of this Gaussian, a^, is varied 
in the range 0.3 - 0.7 MeV, consistent with beta-delayed neutron emission experiments. 
We exhibit the Gamow- Teller strength distribution in Fig. 1 both without and with 
the modulation beyond the particle emission threshold. The results of indicate that 
F^/F for the IAS falls off to 10% at an excitation energy of 1.532 MeV with respect to 
the neutron emission threshold. In the case under consideration the suppression will 
be much larger due to lack of isospin constraints, as noted earlier. For comparison, 
we choose = 0.5 MeV which makes exp(— x^/cr^) = 8.4 xlO~^ at x = 1.532 MeV. 
Modulating the contribution beyond n-emission in this manner, the total capture rate 
is 133.09 (137.46) SNU, an increase by 0.26 (0.30) SNU over the standard value using 
the SSM of ref. |T| (0). The contribution from the different sources are given in Table 
1, where we see that the increase comes almost entirely from ^B neutrinos. We have 
checked that the increase in the capture rate does not depend very sensitively on the 
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parameter cr^: varying from 0.3 to 0.7 MeV changes the total neutrino capture rate 
from 132.99 (137.36) SNU to 133.16 (137.55) SNU using the solar model of ref. 

(i)- 

The work discussed so far assumes that the model space for the distribution of 
the GT strength includes only the four fp shell orbitals, namely I/7/2, 1/5/2, 2^3/2 
and 2pi/2. However, the experimental pick-up reaction data suggest that the l5'9/2 
neutron orbital is about 10% occupied |ll|] in the ground state of ^^Ga and one also 
sees that '^^Ge has a low-lying excited state of 9/2"*" at an excitation of only 0.1984 
MeV [|12|. Also, theoretical estimates using spectral distribution methods and with 



different interaction Hamiltonians in the (2^3/2, I/5/2, 2pi/2, '^99/2) shell model space 
show that the 1(79/2 orbit is about 16 - 20 % occupied |]l3l in the ground state of ''^Ga. 
Inclusion of the lgg/2 orbital in shell model calculations makes the model space too 
large to handle. On the other hand, using our schematic model of the GT strength 
distribution one can extend the calculation to include it rather easily |T^|. Here the 



ground stae of "^^Ga is a pure configuration with the fp neutron orbitals full except 
1/5/2 which has 4 particles instead of 6 and the two extra particles are put in lgQ/2- 
This configuration gives rise to two additional transitions other than the ones listed 
earlier, namely lg<ji2{n) — > lg<ji2{p) and lgg/2{n) —> lgT/2{p). For the GT strength 
distribution this will give rise to seven Gaussians and one half-Gaussian instead of five 



Gaussians and a half-Gaussian considered earlier The centroids of these extra 

strength Gaussians are evaluated in the same manner and they are given the same 
width parameters as others. The f3~ sum rule strength, S13-, increases in this space 
to 221/7 from the earlier value of 3(A^ — Z) - i.e., 27 - giving rise to S13+ of 32/7. 
With this form of GT strength distribution one can again investigate how much the 
neutrino capture rate increases for the Ga detectors if one includes the contribution 
beyond the neutron emission threshold in ''^Ge. Using the same model for the ratio 
of gamma decay width to neutron emission width, one gets the capture rates in this 
fpg model spaces as given in Table 2 for the solar models of ref. 0] and 0. These 
numbers are for the same value of the strength width parameter of 1.55 MeV and 
= 0.5 MeV as in the fp case. We see that the increase in the total capture rate from 
contributions beyond the n-emission threshold is 0.32 SNU for fluxes of ref. [|l| and 
0.36 SNU for the solar model of ref. 0. Varying from 0.3 MeV to 0.7 MeV changes 
the total rate from 133.182 (137.572) SNU to 133.386 (137.806) SNU using the solar 
model of ref. i (§). 

A microscopic calculation of the ratio of the gamma decay width to the total width 
is being attempted. In conclusion, we stress that the contribution to the neutrino 
capture rate of the excited states of "^^Ge beyond the neutron emission threshold are 
not negligible and are actually larger than the Fermi transition contributions usually 
included. Taking into account the contribution beyond the n-emission threshold in 
^^Ge along with the extension of the space to fpg orbits give an increase as large 
as 0.37 SNU over the conventional theoretical calculation done including the fp shell 
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orbits only. 
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Table 1: Neutrino capture rate on ^^Ga without (columns A and C) and with 
(columns B and D) contributions above the neutron emission threshold using 
Bahcall-Ulrich and Bahcall-Pinsonneault solar neutrino fluxes. Here has been 
chosen to be 0.5 MeV. 



Solar Neutrino 
source 


Bahcall - Ulrich SSM 


Bahcall — Pinsonneault SSM 







A 


B 


C 


D 


pp 


67.781 


67.781 


66.727 


66.727 


pep 


5.211 


5.211 


5.211 


5.211 


hep 


0.043 


0.045 


0.043 


0.045 




34.109 


34.109 


37.490 


37.490 




11.335 


11.593 


13.036 


13.333 




4.823 


4.823 


4.823 


4.823 


15^ 


9.433 


9.433 


9.742 


9.742 




0.093 


0.093 


0.093 


0.093 


Total 


132.828 


133.088 


137.165 


137.464 
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Table 2: Same as in Table 1 but for the fpg case. 



Solar Neutrino 
source 


Bahcall - Ulrich SSM [0 


Bahcall - Pinsonneault SSM 


A 


B 


C 


D 


pp 


67.781 


67.781 


66.727 


66.727 


pep 


5.211 


5.211 


5.211 


5.211 


hep 


0.044 


0.047 


0.044 


0.047 


''Be 


34.109 


34.109 


37.490 


37.490 


'B 


11.485 


11.798 


13.208 


13.568 


13iV 


4.823 


4.823 


4.823 


4.823 




9.433 


9.433 


9.742 


9.742 




0.093 


0.093 


0.093 


0.093 


Total 


132.979 


133.295 


137.338 


137.701 
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FIGURE CAPTION 



Fig.l. The Gamow- Teller strength, Bgt, (solid line) as a function of the excitation 
energy of the Ge nucleus. Also shown is the neutron-emission threshold in Ge at 
7.4 MeV. Modulating the distribution beyond the particle emission threshold by a 
Gaussian suppression factor of width 0.7 MeV (0.3 MeV) results in the long-dashed 
(short-dashed) curve. 
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